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By exploiting the relatively long lifetime of the tau lepton, we propose several novel methods 
for searching for new physics at an electron-positron collider. We consider processes that involve 
final states consisting of a tau lepton pair plus two missing particles. The mass and spin of the 
new physics particles can be measured in 3-prong tau decays. The tau polarization, which reflects 
the coupling to new physics, can be measured from the r —>■ irv decay channel using the impact 
parameter distribution of the charged pion. We also discuss the corresponding backgrounds for 
these measurements, the next-to-leading order (NLO) effects, and the implications of finite detector 
resolution. 


I. INTRODUCTION 

Despite the excellent successes of the Standard Model (SM) of particle physics it is clearly incomplete, e.g., it does 
not explain Dark Matter (DM), the gauge hierarchy problem or the current excess of matter over antimatter in the 
universe. Supersymmetry (SUSY) [Tl[2] is one of the leading candidate theories for new physics, i.e., physics theories 
that go Beyond the Standard Model (BSM). If supersymmetric particles (sparticles) have masses below the TeV scale 
then the gauge hierarchy problem in the SM is solved naturally. In supersymmetric extensions of the SM (SSM), 
the gauge couplings are unified at a high energy scale around 0(10^®) GeV [3H5], which is a strong indication of 
the possible existence of an overarching Grand Unified Theory (GUT). With the assumed existence of a conserved 
R-parity in these SSMs the lightest supersymmetric particle (LSP) become stable and so serves as a potential DM 
candidate [5]. 

Gonsiderable effort has gone into searches for evidence of SUSY at the Large Hadron Collider (LHC) at CERN. 
The lack of any evidence for SUSY to date has pushed the potential masses of sparticles to higher scales, which gives 
rise to the so-called little hierarchy problem [71[S] for SUSY. At the 8 TeV LHC with an integrated luminosity of ~ 20 
fb“^ a gluino mass below 1.3 TeV and the first two generation squark masses below 850 GeV have been excluded 
with a 95% Confidence Level (C.L.) in simplified models (HUTU]. Third generation squarks have looser constraints, 
primarily because the proton contains no heavy valence quarks and because the SM tt background is large. The 
exclusion bounds have gone up to < 700 GeV for both stop mm and sbottom mm when the LSP is light. 
Sparticles in the electro-weak sector can still be relatively light m. since the current searches at the LHC can only 
exclude left-handed selectron and smuon masses that are < 300 GeV [min]. There is no current LHC bound for the 
tau slepton (the stau) [18]. The strongest bound for the stau is from Large Electron-Positron Collider (LEP) searches 
is relatively modest and is rrif > 90 GeV [121 . The possibility of a light stau is well motivated in many SUSY model 
frameworks. In Generalized Minimal Supergravity models mm sleptons are much lighter than squarks due to the 
lack of an SU{3)c coupling. Furthermore, the large Yukawa coupling of the tau lepton makes the lighter stau even 
lighter as it causes the stau mass to decrease more rapidly than the first two generation sleptons when evolving down 
from the GUT scale. 

The main reasons for the modest bounds on the stau are its small production cross section and the relatively large 
backgrounds at proton-proton colliders. A future e“'"e“ collider, where the processes are dominated by electroweak 
coupling, will provide a much improved environment for stau searches and would provide an ideal tool for the study 
of the electro-weak sector of SSMs. Indeed, one of the main tasks for the proposed International Linear Collider 
(ILC) |22H24] is to provide excellent track reconstruction with fine momentum and impact parameter resolution in 
order to have a better measurement of Higgs couplings, such as /irr. The impact parameter resolution |25j is expected 
to be (Tr.0 = 5 ^m 0 [I0/(p/GeV) sin^^^ 6] fim where 6 is the angle between the lepton and the beam axis in the 
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lab frame. Note that such precision can only be reached in the rcj) plane, since the nominal size EIIIT] of the beam 
bunch at 1 TeV is 639 nm (width) x 5.7 nm (height) x 300 /rm (length), where the length in the z-direction is much 
larger than the width and height. 

Following the next LHC run it may well be that the stau (and any new particles with similar decay channels) remains 
unobserved even if it lies within the reach of a future e“'"e“ collider. Possible new physics processes with electron 
and muon final states have been studied extensively in the context of an ILC and for the LHC [28H32]. The energy 
distribution of the final state leptons can be used to measure the sleptons’ mass precisely [33], while the threshold 
behavior of the excitation curves and the angular distribution of the observable particles can be used to determinate 
the spin of new physics particles [34] . Final states involving the tau, while more complicated to analyze, are interesting 
and potentially very informative. The end points of the tau jet energy distribution can be used directly to extract 
the masses in the decay chain [35] . while the distribution of the t ^ ttv channel and the E^± /Ep± distribution of 
the p^v channel are able to measure the tau polarization |36H41j . Furthermore, the spin correlation of two tau 


leptons in the H ^ tt channel can be used to determine the CP properties of the Higgs boson |42H45| . Some recent 
studies also show that vertex information for the tau decay can help in the reconstruction of the H —>■ rr |42l |46] 
process and other SM processes m • 

In this study, we will take advantage of the anticipated powerful track reconstruction capabilities of the ILC to 
extract additional information about the final states of the tau decay and use this to assist in event reconstruction and 
the search for new physics, e.g., stau pair production with a subsequent decays of the form r —>■ tx^. The discovery 
prospect of a 1 TeV ILC is analyzed here for a relatively heavy benchmark point of rrif = 300 GeV and m^o = 50 
GeV. However, our method will also be useful at other collision energies. Reducing the collision energy while staying 
in the kinematically accessible region would increase the production rate of our signal. 

The paper is organised as follows. In Sec. |lT| we give a brief overview of tau decay. In Sec. jlllj we present the 
theoretical framework for the mass and spin reconstruction and we propose and discuss some techniques for the 
determination of tau polarization. The parton level results as well as their statistical uncertainties are also discussed. 
A more detailed analysis including detector effects is presented in Sec. IV and finally we discuss our results and give 
our conclusions in Sec. El 


II. TAU DECAYS 

The tau lepton has complicated decay branching fractions. Its dominant decay channels and their branching 
ratios |48| are given in Table |T] The 27r and 37r decay modes are dominated by the p and ai meson resonance 
respectively. Based on the number of charged particles in the final state, the hadronic decay modes are classihed into 
two classes referred to as 1-prong tau decay and 3-prong tau decay. In this work, only those tt and p/ai resonance 
channels which are shown explicitly in the table will be considered. Note that those contributions comprise around 
85% of total hadronic tau decays. 


Leptonic decay 


1-pror 

p^(—7r*7r°)j/ 

ig decay 

ai(— 7r^7r°7r°)i/ 

others 

3-prong decc 

fli (—>■ 7r*7r^7r^)i/ 

'■y 

others 

35.2% 

10.8% 

25.5% 

9.3% 

6.4% 

9.0 % 

3.8% 


TABLE I: Tau decay channels and their branching ratios. 


Since we cannot measure the momentum of the neutrino, we will not be able to fully reconstruct the tau momentum 
in the general case. However, the tau provides more information about potential new physics interactions than the 
first two generation of leptons through the kinematics of its decay products. The polarization of the tau is correlated 
with the energy ratio between the energy of the visible decay products and the energy of the tau lepton as shown in 
Fig.[T] The correlation between the energy ratio and the tau polarization depends on the decay products, where we 
see that the strength of the correlation decreases as we go from tt to p to ai. However, it was shown in Ref. [IHI that 
by considering more complicated multi-dimensional distributions, stronger correlations with the tau polarization for 
the p and Oi channels can be found. 

References (33] 1351130] were devoted to the reconstruction of the mass and the coupling of the stau based on the 
T ^ pv channel. In those references the mass reconstruction was seen to suffer from limited statistics, especially 
when extracting the endpoint of the p meson energy distribution. In the present work we shall try to use another 
remarkable feature of the tau, its relatively long lifetime, to obtain a more precise reconstruction of new physics masses 
and couplings and of the spin of any new physics particle. The tau lepton has a mean lifetime of r = 2.9 x 10“^^ s 
and a mass of 1.777 GeV. The probability of a tau traveling further than length L is P{j3) = with Lq = 
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FIG. 1: Energy ratio between visible and total final states for tl (left panel) and tr (right panel) decay. The 

distributions are calculated by TauDecay [50] . 


"fCT = 




= —CT . So, a typical r lepton with energy ^ 0(100) GeV will travel a distance I ^ 5 mm before 


its decay. If the tan goes through a 3-prong decay, the reconstructed secondary vertex will be displaced from the 
primary vertex of the signal process. From the position of the displaced vertex, we can obtain the direction of the t 
lepton. Ref. |46j showed that with the information obtained from the displaced vertex there is an improvement for 
both Higgs detection and the reconstruction of the Higgs. A similar technique is used in this work to reconstruct the 
more complicated process of stau production with a subsequent decay into a tau lepton and a neutralino. We will 
show that with knowledge of the tau direction, the corresponding tau energy can be reconstructed. This can be used 
to learn about the masses of new particles. If the directions of both tau particles in the final state are known, the 
whole system can be resolved. As a result, the angular distribution of the stau can be used to determinate its spin. 
We also propose a new method for studying the tau polarization through the impact parameter distribution of tt in 
the T ^ TTv channel. This method benefits from a higher sensitivity to tau polarization m and is potentially more 
accessible experimentally, since it does not require the reconstruction of a p candidate. 


III. NEW PHYSICS RECONSTRUCTION 


In this section we consider the following two benchmark processes in SUSY models that can gives rise to a 2r -|- 
signature. They differ in the spin of the new particles: 


e+e t+t + X°X° ^ jij 2 + $t (1) 

e+e" xtxi t+t~ + vv ^ jij 2 + $t , (2) 


where ji _2 stand for the visible components of the r decays and /v are invisible for the detector. The chargino 
benchmark is studied for the purposes of making a comparison in spin reconstruction. It will not be considered as a 
background for f searches, since our study will be restricted to the consideration of only one new physics particle at 
a time. Both processes have contributions from s-channel 7 /^ mediation. The second process may have additional 
contributions from a t-channel Vg mediation, which is suppressed in the case of using a right handed electron beam 
or a heavy sneutrino. Since a right handed electron beam is considered in our study in order to suppress the WW 
background, only the s-channel contribution needs to be considered. The same Feynman diagrams can show up in 
many other new physics frameworks as well, e.g., in models with universal extra space dimensions [521153] . 

In the process represented by Eq.Q, the polarization of the produced r is determined by the degree of mixing of 
the stau sector and the neutralino sector. The bino-like neutralino gives for the tau polarization 


B(rnX°) - B(tlX°) 
BiTRX°) + B{tlX°) 
( 4 sin^ Q-t — cos^ O-r 


4 sin^ Qr + cos^ Qr 




( 4 ) 
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where Or is the mixing angle for the stau sector. The polarization Pr in this case can lie in the range [—1,1] when 
varying Or- Only the tau polarization will be of concern in our study and so the model with a bino-like is 
representative of a large class of models with arbitrary neutralino mixing. We will use this model framework to study 
different tau polarization effects. The Monte Carlo events are generated using MadGraphb [M] with the TauDecay [50] 
package used to perform simulations of the decaying tau lepton. 


A. Knowledge of the tau direction from the displaced vertex 


For the two processes given above, there are two visible jets and four missing particles in the final states. At first 
sight it might seem impossible to fully reconstruct the decay system. However, as mentioned in the introduction, 
additional information about the tau lepton can be extracted precisely because of the relatively long tau lifetime and 
the consequent displaced tau decay vertex. 

With the measurement of the r direction from its displaced vertex, we can write the t four-momentum as 

p'f = {Er, \pt I sin 0 cos (j), \pr I sin Osincf), \pr \ cos 9} (5) 


and we can measure the four-momentum for the associated jet 




( 6 ) 


where Er = \/|_Pr p + Using the fact that the neutrino is almost massless, = {pr — Pj)"^ = 0, we can then 
solve the kinematics of the system up to a two fold ambiguity 


-2CD ± - 4(£)2 _ AE^){C^ - AEjml) 

2{D^-AE'jy^ 


(7) 


where 


C = ml +m'j (8) 

D = 2pl sin 9 cos (j) + 2p| sin 0 sin (/) -|- cos 0 . (9) 


In summary, from the jet measurements and the direction of the displaced r vertex we can reconstruct the r four- 
momentum, up to a two-fold ambiguity. 

Prior to the r decays our processes of interest contain two t’s and two missing neutral particles as indicated in Eqs. 
Q and ([^. We use the following notation in the reconstruction of the full process 


Pr^ = {e^\pI\pI\pI^) , 
Pu, = {E^\p2\Py\pT) > 
Pr2 = {E''^,P?,Py',P?) , 
Pn2 = {E”-\Px\Py\PT) ■ 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Up to the two-fold ambiguity for each of the two r four-momentum there are eight remaining degrees of freedom. 
However, we have the following eight constraints from four-momentum conservation and on-shell mass conditions in 
the centre of momentum frame. 


E^l ^ , (14) 

p -b p -b p -f p = 0 , (15) 

{E-^f - iP^r - {P7 f - {p7f = ml , (16) 

{E-^f - {p7? - {P7f - (pT? = ml , (17) 

(ijn + )2 _ (px, _ pm )2 _ (px, + pm )2 _ (pm + pm )2 ^ ^2 
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By solving this system of equations we find 


where 


^ i = 1,2 

p "2 ^ _^ri _ ^T2 _ -ni 

=D, + D2p;^ + DspT 
p^^ =F,+f2p:^ 

^ -{G 1 G 2 + F 1 F 2 ) 

P- Gl + F| + 1 

v/(GiG2 + F,F2r - (Gj + Fj + l)(m2 +Gl + Ff 

Gl+Fi + 1 


(£;ni) 2 ) 


Cl = 

Cl p'ri+'r2 j,Ti+T2 

n _ 1 n _ _ r\ _ _ t’z 

1 9„Ti+r2 ’ -^2 ^^1+^2 ’ 3 „ti+T2 > 

Px Px 

El = 2E^^E^^ - 2pl^Di + {rrE^f + {nFf - {m^ f, 
E 

^ 2{pI^D2+pI^)’ 

p p7D3+p7 

^ pV-D 2 +py-' 

Gi=Di+D 2 Fi, 

C2 = D2E2 + -D3 . 


(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 

(24) 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


There is a two fold ambiguity also for this system of equations. 

By reconstructing the kinematics in this manner, we will have on overall eight-fold ambiguity made up as 
2(from Ti) X 2(from T 2 ) x 2(from final reconstruction) = 8. We will next show that the new particles mass and 
spin can be reconstructed by using this method. 


1. Reconstructing new particle masses 


For the processes described in Eqs. 0 and 0 the 
end points given by 


energy distribution is bounded from above and below with its 


El 


max( min 


(min) 


E;±piPf 


(32) 


where 

mio + ml ,_ 

-, Pr = V(.E*y -ml, /3f = (1 - 4 mj/s)'^G _ ( 33 ) 

By studying the energy distribution of the tau lepton, we will be able to reconstruct the mother particle (f or x^) 
mass and missing particle {x^ or v) mass. 

With the four momentum of the visible sector from the r decay and the direction of the r momentum, the energy 
of the r can be reconstructed up to a two-fold ambiguity. Up to small changes due to next-to-leading order (NLO) 
corrections and detector effects, we know from two body decay kinematics [33] that the distribution of the “true” 
tau energy solution is flat over the entire allowed energy range, Emin < E < Emax, and independent of the tau 
polarization. In general we expect the energy distribution of the r decay products to depends on the tau polarization, 
which may lead to different distributions of the false Er solution for different tau polarization. By considering the 
r pair production process we show the distribution of Efaise and Eaii for three different tau decay channels in Fig[^ 
where we have taken = 300 GeV, m^o = 50 GeV and y/s = 1 TeV. 

From the figure, we can conclude that for the t ^ ttv and the t ^ pv decay channels, the false tau energy solution 
distributions are strongly dependent on the tau polarization. As a result, it will be difficult to extract the true tau 
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FIG. 2: Reconstructed r energy for different chiralities and decay topologies: r —>■ tti/ (upper left), t —>■ pi/ (upper 
right) and r — > aiv (lower left). The lower right panel gives the full tau energy distribution for mf = 300 GeV and 
myo = 50 GeV. The arrowed lines show the reconstructed endpoints (i?max/min) for Etme distributions. 


energy distribution without knowing the tau polarization. The situation is different for the t — >■ aiv, which gives 
the dominant contribution to 3-prong tau decay. For this channel, different tau polarizations give almost identical 
false tau energy distributions. Moreover, the tau direction in this case can be reconstructed from the location of the 
secondary vertex. 

For illustration purposes, a sample of 3000 three-prong tau decays is used to study the mass reconstruction precision. 
The superposed distribution of Fltrue and iSfaise are shown in the lower right panel of Fig. with statistical uncertainty 
built into the distribution. The full (true-|-false) distribution is comprised of a rectangular i^true distribution and a 
continuous Flfaise distribution. We use following algorithm to locate the endpoints of the i?true distribution from the 
full distribution: 


• The location of the falling edge is estimated to be where (iVj_i — is maximized, where W is number 

of events in ith bin; 

• The height of the rising edge which is given by the height of the rectangular i?true distribution, can be estimated 
by /i = W_i - W+i; 

• Allowing some fluctuations, the largest j such that Nj — Nj-i > 0.8h can be used to locate the rising edge. Note 
only two adjacent bins are used here because of the better statistics due to the shape of the false distribution; 

• Improved estimates of the location of the edges are given by by Amax = + f Amin = 

fo-i-fo-i , where Ei is the center value of the tth bin and S is the bin size; 


2 
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The corresponding uncertainties are = S/y/S and 


hs/VT2)^+ Y. 


,dE„ 


1,2,2+1 


dNk 




(34) 


The reconstructed central values of i?niax and ^'mm 


are shown by arrowed lines in the lower right panel of Fig. 
for the left-handed tau (i?min = 50.0 ± 3.5 GeV, ifmax = 434.4 ± 3.2 GeV) and the right-handed tau (ifmin = 50 ± 3.5 
GeV, i^ax = 438.5 ± 3.5 GeV) respectively. The new particles masses can then be reconstructed by using Eqs. 


and (|33|). For left-handed tau we obtain = 304.2t®g GeV, m^o = 53.8l(g GeV and for right-handed tau we find 


- 1-6 


mf = 303.llg GeV, m^o = AG.O'^l GeV, which compare well with our input values of rrif = 300 GeV, myo = 50 GeV. 


2. Reconstructing new particle spin 


For the s-channel y/Z mediated scalar/fermion production at an e“''e“ collider, the spin of the new particles can 
be related to threshold excitation and angular distribution |34| . Because of our relatively heavy benchmark point, 
accumulating a sufhcient number of events at different collision energies would require a large amount of accelerator 
operation time. Therefore, we will focus exclusively on the angular distribution as a spin discriminator in this work. 
The polar angle distribution for the scalar pair production in the s-channel in the Gentre of Mass Frame (CMF) is 


1 da 


a d cos 0n 


[e+e 


f f ] oc sin^ 0 - = 1 — cos^ 9f 


(35) 


whereas for fermion pair production in the s-channel we have 


1 da 


a d cos 9y^± 


[e+e 


oc (s - 4 to|± ) cos^ 9^± -I- (s -I- 4m|±) . 


(36) 


From these equations, we can see that a scalar particle tends to be produced in the central region, whereas the 
fermion is more likely to be produced in the forward/backward regions. It can be shown that the t-channel v mediated 
process for fermion final states will lead to an asymmetric distribution for cos 9^± when only one particular charge is 
considered [^. For each electric charge of the final state, either the forward or the backward direction is favored. In 
particular, when the sneutrino mediator is light, the final state is much more concentrated in the forward/backward 
regions than the s-channel process, which helps make the spin characteristic even more distinguishable. However, the 
t-channel process is highly suppressed for a right-hand polarized electron beam and may even be non-existent in some 
new physics scenarios. For this reason, we do not consider this subprocess in our discussions. 

Using Eqs. (35) and (36), we can reconstruct the polar angle distribution for t/x^ in order to extract spin infor¬ 
mation about the new particles. As we have noted before, our method typically produces an 8 -fold ambiguity for our 
system with only one of them being physical. The experimentally accessible variable is the 8 -fold superposition of 
all solutions. In practice, the distribution of the seven false solutions will inherit information about the true physical 
solutions and will be different for scalar and fermion particles. We will restrict our considerations to the case where 
there is only one new physics process at a time. For example, we do not consider here the possibility of a superposition 
of the effects of both scalar and fermion new physics particles. 

The polar angular distribution of the sum of the 7 false solutions (denoted as False) and all 8 solutions (denoted 
as True-|-False) for f and are shown in Fig. In the figure we include both charges of the produced particles 
(f, x^), which will lead to a symmetric distribution for the t-channel process. From the figure we find for all cases 
that the corresponding false solutions are distributed similarly but are flatter than the true solution. The t-channel 
contribution is shown as well, where the sneutrino mass is taken as 100 GeV. It is featured by a more remarkable 
concentration at forward and backward region than s-channel fermionic final state process. 

In order to measure the difference between particles spins, we can define a spin measure on the superposed distri¬ 
bution of all eight solutions (corresponding to True-I-False): 


A., = 


N{\ cos(0)| < 0.5) — N{\ cos(0)| > 0.5) 
N{\ cos(d)| < 0.5) -I- N(\ cos(6*)| > 0.5) 


(37) 


We see from Fig. that Ag is positive for scalars and negative for fermions. As discussed above, we will only consider 
the s-channel process here. Using a very large number of simulated events we construct the red lines in Fig to high 
precision. We can then extract Ag. We refer to this as the theoretical prediction for the spin measure and denote it 
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FIG. 3: The cos 9 t distribution for f (left panel) and chargino (right panel) pair production at a 1 TeV e“''e“ 
collider. The red line and green shadowed region in both plots correspond respectively to the True+False and the 
False solution distributions for the s-channel process. The blue line and yellow shadowed region in the right plot are 
respectively the True+False and the False solution distributions for the t-channel fermion production. 


as We find ^fermion = —0.065 and = 0-2- The corresponding statistical uncertainty [56] can be estimated 

by 





N '■ 


(38) 


where iV< and 7V> stands for iV(|cos(0)| < 0.5) and 7V(|cos(0)| > 0.5) respectively. The statistical uncertainty of 
^<{>) Because there is a correlation between those 8 -fold solutions, the total number of events 

N is used to give a conservative estimate for the statistical uncertainty rather than 8 x iV. As a result, without taking 
into account any detector effects, an event sample number of 


N 


9 X 


1-KLlar)^ 


_ 4 th 
'^^scalar ^fer 


fermion 


120 


(39) 


is sufficient to distinguish the stau spin at the ia level. 


B. Tau polarization from impact parameter 


Tau polarization can be inferred in all decay channels in different ways. The simplest case is r —>■ ttz/ with the left 
handed interaction m 

C = V2Cf^^PLi'rdfj.T^ + h.c. . (40) 

In the tau rest frame, the relation between tau polarization(PT-) and the pion angular distribution(cos 0) is given by, 


- = -(I 

F d cos 6 2 ^ 


Pr cos 0). 


After boosting the whole system to the Laboratory frame (LAB), we will have 

+ ml + j3^p^ cos 9 , 

Er = jrrir , 


(41) 


(42) 

(43) 


where 7 is the boost and where = (m^ — rn^)/2mT and 9 are the magnitude of the pion three-momentum and its 
polar angle in the tau decay rest frame respectively. The energy fraction of the pion in the LAB frame x = E-j^jEr is 
linearly related to the pion polar angle in the tau rest frame for a highly relativistic tau (/3 —^ 1 ), since we have 


Ett 

Et 


cos 9 /3-^l 


= (i + :Lf)h+/3(i-wf)- 


(1 


mi 2 mi 


X = 


2 


2 


(44) 
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So, from Eq. (41) we find in the Lab frame after a change of variable from cos 9 to x 


1 dr 1 2x — {1 + {rriTr / nir)'^) 

r da; 1 — (m-^IvcirY ^ (1 — [rriT^lmrYY 

- 1 + P^(2a;- 1), 


(45) 

(46) 


which can be used for studying tau polarization [36]. 

However, in the 7r-channel, the energy of r lepton cannot be easily reconstructed for most processes of interest, since 
there is usually more than one missing particle in the final state. Exceptions arise in some very special cases, e.g., (i) 
for a few TeV tau lepton whose energy can be measured by its track curvature directly, (ii) for single tau production 
with a hadronic decay, and (iii) for tau pair production at electron-positron colliders. Those exceptions are beyond the 
scope of our current study. The tt energy spectra can also be used directly to measure the tau polarization for known 
new physics process [311 [ST] HD] . However, the tt energy distribution depends on the process under consideration and 
on the masses of new physics particles. 

In most studies [SHlIMlIll] at colliders the r —)■ p(—7r+7r°)i/ channel is chosen for measuring the tau polarization due 
to the fact that a decays mostly to a longitudinally polarised p meson while a decays mostly to a transversely 
polarised p. The energy ratio z = E^± /Ep can be used to measure the polarization of the p meson, which in turn 
gives information about the tau polarization. The distribution of z with respect to the tau polarization is given in 
Ref. [36| . Their equations show that the z distribution is related to the tau polarization in a complicated way. Studies 
in Ref. [51| have given the general relative polarization sensitivity for each tau decay channel. They found that the 
T ^ pu channel does not perform as well as the r —>■ ttiz channel, while the three prong decay and leptonic decay 
channels only have 1/3 the sensitivity of the r —> ttz/ channel. Another difficulty arises for the p channel in practice. 
One needs to have an efficient tagger on the p meson in order to recognise the decay topology since both mis-tagged 
tau jets and other 1-prong tau decay modes have to be excluded first. 

The proposed ILC has the advantage that it can resolve the track precisely. As mentioned in the introduction, an 
impact parameter resolution as small as 5 /rm can be reached in the rcj) plane. This will be very helpful if we try to 
study the tau polarization using the r —>■ -ku channel. Eor those processes of interest, the tau lepton energy has a flat 
distribution between Emin and Emax- Moreover, different tau polarizations give rise to very different distributions of 
the pion polar angle in the tau rest frame. This means that in the laboratory frame we will have a different impact 
parameter distribution for a different tau polarization, which in turn is related to the angle difference between the 
tau and pion momenta, 9^_t. The impact parameter is given by 


d = LsuiOtt-t , 


(47) 


where L is tau decay length. Therefore the impact parameter is dependent on the tau polarization. 

A difficulty arises in a realistic detector. Because of the long lifetime of the tau, there are no long-lived hard 
tracks associated with the interaction point (IP). In addition, the interaction region is relatively wide along the z- 
axis Therefore we can not have good resolution for the position of IP. As a result, the impact parameter becomes 
experimentally unreconstructable. Thanks to the narrow beam size along the vertical and horizontal directions, the 
position of IP can still be measured precisely in the rcj) plane. Then, we can define the geometrically signed impact 
parameter m which is the distance of closest approach of the extrapolated track to the assumed production point, 
the centre of the interaction region, in the rcj) plane, 

dcip = L sin 9t sin(^ — (/r) • (48) 


In the following, we shall study the sensitivities of variables to the tau polarization. Eor the t —>■ iri) channel, we 
have shown in Eq. (44) that the energy ratio (x = Et^/E-t) in the LAB frame has the same sensitivity as the polar angle 


6, since they are linear related with each other in the collinear (i.e., relativistic) limit. It should be noted that neither 
X nor 9 are reconstructable experimentally for the stau pair production process, except in special circumstances as we 
noted earlier. However, we will use x extracted from the Monte Carlo data as “best-case reference” for the purposes of 
comparison with the impact parameter dcip- We first study the x = E^/E,- distribution for different tau polarizations, 
focusing on three benchmark points, Pr = I.O, Pr = 0.2 and Pr = —1.0. We use these to demonstrate the sensitivity 
to right-handed, highly mixed and left-handed polarized tau, respectively. The variation of the polarization by ^ 0.1 
around these three benchmark points is also studied for comparison. The x distributions are given in the left panel of 
Fig. for all those 7 cases. In addition, the Monte Carlo data with limited statistics (3000 tau decays) for those three 
benchmark points are superimposed on the plot. The error bars show the corresponding statistical uncertainties. To 
have an intuitive comparison with other variables we define the following measure for tau polarization. 


Artt = 


N. 


x<0.5 


-N. 


x>0.5 


K 


x<0.5 


+ w 


a:>0.5 


(49) 







10 


which is sensitive to the shape of the dependence on x. As for the t ^ pv channel, we show the energy ratio of 
z = E^± /Ep in the right panel of Fig. where a similar analysis to that used for ttv channel has been used. From 
the figure, we can define another measure of the polarization sensitivity for z, 


Arp = 


-^|z-0.5|<0.25 ~ A^|z-0.5|>0.25 
-^|2-0.5|<0.25 + A^|z-0.5|>0.25 


(50) 



0.4 0.5 



E,±/Epj 


FIG. 4: Theoretical and MC truth distributions by using discriminate variables constructed from t ^ ttr (left 
panel) and t —>■ pv (right panel). The MC truth distribution is using 3000 tau decays with statistical uncertainties 

of for each bin i . 


Finally, we propose to use the information from the impact parameter of the pion track in the r —> tti^ channel to 
study the tau polarization. This variable is more easily accessible experimentally than z = Ej^±/Ep. The decay length 
of the tau is generated with the probability distribution P = . Hence, we can calculate the geometrically 


signed impact parameter by using Eq. (48). The distribution of the geometrically signed impact parameter ^gip is 
shown in Fig. where the curves for all polarizations of taus appear to be intersecting at around dcip ~ 41/im. So 


we can define 


. -^|dGip|<41 - A^|dGip|>41 

>lriTP =- — 


-^|dGip|<41 + A^|dGip|>41 


(51) 


as a measure for studying the tau polarization for this variable, where fV|dGip|<4i i® number of events with 
MgipI < 41/xm. Note that the modulus is necessary since the definition of c^gip allows for it to be negative. A 
remarkable feature of dcip as shown in the left panel of Fig. [^is its insensitivity to the energy of the tau lepton. A 
simple explanation of this is that for an energetic tau, the decay length extension due to the additional boost along 
the tau momentum direction is canceled by the narrowed angle difference between pion and tau. As a result, for tau 
energy larger than 0(10) GeV, the distribution of cIgip is strongly dependent on the tau polarization and very weakly 
dependent on the tau energy. For tau energy lower than ^ 5 GeV, the energy dependence of c^gip increase. In this 
region, the tracks from tau decay tend to have energy smaller than 1 GeV which can not be reconstructed effectively 
at ILG and tau identification may also suffer from heavy contamination. However, as we can see from Fig. there 
are few events with tau energy in this region. 

Having defined measures for all those variables, we are able to study their relative sensitivities to tau polarization. 
Firstly, we will calculate the corresponding measure for each polarization with a large number of simulated events, 
where the statistical uncertainty is very small. We refer this value as the “theoretical prediction”. We also consider 
measures for a much smaller Monte Garlo data set (containing 3000 tau decays) which is referred as the “experimental 
measurement”. Then, the corresponding statistical uncertainty for the smaller data set can be calculated by using 
Eq. (38). The results are presented in Fig. 

From Fig. we find that the experimental measurements agree with the theoretical predictions within uncertainties 
in all cases. For a given Pr, the sensitivity of each measure is proportional to the ratio between the slope of the asym¬ 
metry from the theoretical prediction and the uncertainty of experimental measurement at Pt. Since the differences 
in the slopes of asymmetries are larger than the differences in uncertainties, the sensitivities can be simply estimated 
as the magnitude of slopes. Then the asymmetry Art^ of the x distribution is giving the most sensitive probe as 
expected, but recall that it is not experimentally reconstructable. The asymmetry Agip of the impact parameter is 
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FIG. 5: The polarization dependence of doip for given different energies (left panel) and for a concrete process with 

ruf = 300 GeV and m^o = 50 GeV (right panel). 


Asymmetries of tau polarization measures 



FIG. 6: The three asymmetries as function of the tau polarization Pt. The lines and points correspond to 
“theoretical prediction” and “experimental measurement” as defined in the text. In case of the tau energy has 

been cheated from MC truth. 


less efficient than asymmetry Afip. On the other hand, in practice p meson tagging is generally used to identify the 
T ^ pv signal and to suppress backgrounds. However, this will weaken the sensitivity of A^p. A tau decaying into a 
single charged tt has a higher tagging efficiency than a tau decaying into a p meson. All those considerations suggest 
that the impact parameter of r —> ttiz decay as a very promising measure for the study of tau polarization. 


IV. SIMULATION AND RESULTS 


In the previous section, we have reconstructed the new particle mass and spin by using extra information from the 
secondary decay of the tau lepton. We also propose a new method using the impact parameter to measure the tau 
polarization which can provide information about the new physics coupling. 

In a realistic experiment, we have to recognize the features of our signal processes so that the corresponding 
backgrounds can be specified. A relatively pure signal sample can be collected after the backgrounds are subtracted. 
Moreover, a detector can only have finite resolution on energy and direction measurement, which will lead to further 


complications. As has been seen in Sec III we will encounter quadratic equations when we try to solve the system. 
The effects of finite resolution may lead to complex solutions for those quadratic equations, which means that the 
reconstruction has failed. In what follows in this section, the background subtraction, NLO effects and detector 
smearing effects will be discussed. 
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A. Backgrounds and parton shower effects 

Our signal processes are featured by two tau leptons and relatively large missing energy (> 50 GeV) in the final 
state. The channels of interest are t/jT; + and ThTh + $t- Our study requires at least one hadronically decaying 
tau. The pure leptonic decay mode will not be of interest to us. 

There are many SM processes that could potentially contribute backgrounds in our analysis. However, some of 
them can already be understood to be small with appropriate treatment. SM processes with quark final states can be 
effectively suppressed by requiring low track multiplicity on jets [ID]. It is known that the 77 background is negligible 
when the tau lepton becomes energetic |39j . In addition, a moderate cut on missing transverse momentum frp > 15 
GeV can greatly reduce the eell backgrounds [DSj . 

As a result, the backgrounds of our signal at an e+e“ collider are dominated by WW tutu and ZZ ttuu. 
The Z-boson mediated s-channel and electron neutrino mediated t-channel contribution to WW production and 
the electron mediated t-channel contribution to ZZ production can be highly suppressed by using a right-handed 
polarised electron beam in the collision. The corresponding production cross section for beam polarizations P^,- = 0.8 
and Pf,+ = —0.2 is shown in the left panel of Fig. Some representative cross sections for processes of interest are 
also shown. Those LO cross sections are calculated using MadGraph5 |S1] with the gauge bosons decaying into taus 
and neutrinos. The tau lepton in the final state is required to have pr > 10 GeV and | cos0| < 0.98. As can be seen 
from the figure, the WW backgrounds are smaller than the signal processes of interest by more than one order of 
magnitude in most cases, while the ZZ background is always negligible. 




FIG. 7: Cross sections at a 1 TeV e“'"e“ collider for an electron beam with 0.8 polarization and a positron beam 
with -0.2 polarization. The left panel shows the production cross section for signals and backgrounds. Three 
different stau (chargino) mixings are considered. The right panel gives the tji pair production cross section with 
varying cuts on the transverse momentum of the initial state radiation (ISR) photon. 

As for NLO effects, those real corrections (radiation of an extra photon) can effect our reconstruction methods 
proposed in Sec. m while virtual corrections will only lead to an overall normalization. For the Pr measurement 
the effects of initial state radiation (ISR) lead to a reduced y/s and can be described by a correction factor [ID]. In 
the following, we will show the smallness of such effects in much simpler way. At 1 TeV ILG, the typical f energy is 
^ 500 GeV. The radiation of a photon with energy ^ 0(1) GeV can only affect the stau energy by < 1%. Moreover, 
the production cross sections of a purely right-handed stau {tr) pair with varying cuts on the transverse momentum 
of an ISR photon are presented in the right panel of Fig. From the figure, we observe that there are no more than 
10% of events with an ISR photon energy larger than ~ 0(1) GeV. This leads us to conclude that photon radiation 
effects are negligible comparing to detector resolution that will be discussed later. 

From the discussion above we observe that the backgrounds and the NLO corrections are around one order of 
magnetite smaller than the LO signal process. So, we can safely discuss our method at LO order without considering 
any of those effects. Moreover, other features of our signal may help to suppress the backgrounds even further, while 
keeping our signal intact, e.g. > 50 GeV. However, the precise value of such a kinematical cut depends on the 
model being considered and is beyond the scope of our current work. 
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B. Detector effects 

With the arguments that known backgrounds and NLO effects can be ignored we are then ready to discuss the 
effects of finite detector resolution ^. Let us first note that the tau momentum direction is determined by the direction 
from the IP to the reconstructed secondary vertex. Its precision is mainly limited by the beam bunch length along 
the z axis The position of the IP can be calculated using a technique similar to that proposed in Ref. |58) . 

Firstly, we calculate the probability distribution of the separation angle between the r and total visible final states 
for given r energies, by using the process e“'"e“ Z ^ with fixed center of mass energy. Note that only the 

3-prong tau decay channel is of concern for mass and spin reconstruction in this work. As we can infer from Fig. 
the Ej/Er distribution and thus the separation angle distribution for 3-prong tau decay is not sensitive to the tau 
polarization. So, the probability distribution Pr{Er, A9{T,j)) deduced from above process can be used for all other 
processes with arbitrary tau polarization. Moreover, the e“'"e“ Z ^ process has another two advantages: 

the energy of the r from this process is simply equal to the beam energy and it has a large production rate at the ILC 
which enables precise construction of PriE-r, A0(r, j)). The separation angle distributions for different r energies are 
given in Fig. 


Z ->xi 



FIG. 8: Distribution of the separation angle between r and its visible final states for different tau energies. 

Secondly, in our case, because of the narrow beam size of the ILC, we only have one degree of freedom, which is 
the IP position along z axis, i.e., we denote the reconstructed IP as being at (0,0, z) with the real IP defined to be 
at (0,0,0). For a given z, with the information available from the visible r decay products, we can solve for both the 
Er and the A9{T,j). Note that there is a two-fold ambiguity of Er and the one that provides larger Pr{Er,A9{T,j)) 
is chosen. Two different likelihood functions are defined, Pi(z) = Pr{Er, A9{z)) for a single 3-prong tau in mass 
reconstruction and P 2 {z) = Pt^{Et^, A9{z)) x Pt^{Et^, A9{z)) for two 3-prong taus in spin reconstruction. The z^ax 
which maximize the Pi(z) or P 2 {z) is used as the position of the reconstructed IP for each event. 

In order to estimate the uncertainties of the above reconstruction, we select our f pair events with the r lepton 
energy chosen in two ranges as an illustration, i.e., [49.5,50.5] GeV and [199.5,200.5] GeV. The results are given 
in Fig. where the ag is the angular separation between fo = (Lsin 0,-cos ((),-, L sin0,-sin (/),-, L cos 0) and = 
{Lsm9r cos (j)r,L sin 9r sin (j)r,L cos 9 — Zmax)- We can conclude from the figure, for the single likelihood Pi{z), that 
the reconstructed tau direction is centered on the true tau direction, with l-cr deviation approximated by This 
result is also consistent with Fig. [^ where the angular separation distribution of tau decay indeed shows the width 
around As for the double likelihood P 2 {z), the combinational effect can improve the angular resolution to jp-- 

All other detector effects considered in this work are listed as follows 

• The tracks used in following analysis are required to have Atrack > 1 GeV and lie in | cos6*| < 0.98; 

• Energy smearing for hadronic tau jets is taken as ^ = 3%; 


^ A comprehensive discussion of detector effects for stau search at ILC can be found in Ref. HSl 

^ When there is a small crossing angle(14 mrad) between two beams at ILC, a narrow beam bunch might allow us to reconstruct the IP 
with precision of 0(10) fim along 2 direction. 
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E, =50 GeV E, =200 GeV 




FIG. 9: Angular resolution for three prong decaying tau with Er = 50 GeV (left panel) and Er = 200 GeV (right 
panel). The black curve is for the single likelihood Pi{z) and the red curve is for double likelihood P 2 {z)- 


• Impact parameter resolution in ri^-plane is 5/im0 pjgev] sin^/^ e ’ 

• The tau identification efficiency is assumed to be 0.7. 

From Fig. we determine that the production cross section for the signal process should be at least ~ 10 fb. We 
will take the worst case and work at a 1 TeV e+e“ collider with 1000 fb“^ integrated luminosity to illustrate the new 
physics reconstruction with detector resolution effects included. 


1. Mass reconstruction 

In order to reconstruct the masses of the new physics particles, we will need at least one tau that goes through 
3-prong decay. As a result, the corresponding number of 3-prong tau decays after tau tagging and branching ratio 
suppression is 


~ 10/6 X 1000/6-1 X 0.7^ x (13% x 2) - 1200 . (52) 

The tau energy for mass reconstruction is given by Eq. Q. A finite resolution will lead to negative square root 
in some cases, which means the reconstruction has failed and so the event should be dropped. With the resolution 
parameters given above we find around 1/3 of the total events fail the reconstruction. So we are left with ~ 900 
three-prong tau decays for mass reconstruction. The superposed At^ue/faise distribution for those events are displayed 
in Fig. 

Due to the smaller number of events and the smearing effects of detector, the rising and falling edges are spread 
over more bins than the ones discussed in Sec. IIII A 1 1 Therefore information from more bins needs to be used to 
locate the edge precisely. The modified algorithm is then described as follows: 

• Five bins are used to locate the falling edge, i.e., find the bin where {Ni _2 + A^i-i — A^i+i — A^i+ 2 )/v^ is 
maximized; 

• The height of the Etme distribution is estimated hy h = A /-2 — Ni+ 2 ', 

• Note that the height of Etme will be slightly underestimated because of the smearing effects. So we look for the 
largest j such that A(j+i — JVj-i > h (rather than 0.86, used earlier) to locate the rising edge; 

• Improved estimates of the locations of the edges are given by Amin = A,- -|- f and Amax = 

p S 2Ni+Ni+i+Ni.i-2Ni+2-2Ni.2 . ' ^ 

* 2 Ni-2-Ni+2 ’ 
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• The corresponding uncertainties are 


{S/VU)^+ Y. 

V k=j-i,j,j+i 

SEma. = , /(5/VM)2+ ^ (54) 

V k^i- 2 ,i- 14 ,i+lA +2 ^ 

This gives £^min = 49.2 zb 8.0 GeV, ^max = 444.2 zb 6.8 GeV for the left-handed tau and T^min = 63.4 zb 6.3 GeV, 
^max = 441.1 zb 8.2 GeV for the right-handed tau. This corresponds to the reconstructed rrif = 299.6^22 GeV and 
= 34.5^315 GeV for left-handed tau and to nif = 320.ll^^ GeV and = 4 O.OI 40 GeV for right-handed tau. 
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FIG. 10: The full tau energy distribution for rrif = 300 GeV and m^o = 50 GeV with finite detector resolution. The 
arrowed lines show the reconstructed endpoints (F’max/min) for F^tme distributions. 


T mass reconstruction 
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8. Spin reconstruction 


To reconstruct the new particle spin, both of the taus in the final state are required to go through 3-prong decays. 
So the corresponding accessible number of stau pair events is 


N, = 10/6 X 1000/6-^ X Q.T x (13%)^ ~ 83 


(55) 


We can see from the derivation in Sec. Ill there will be a total of 3 quadratic equations in the reconstruction, which 
means more severe event loss can be expected from complex solutions. On the other hand as a partial compensation, 
and as discussed above, a better IP position resolution can be achieved when we have two 3-prong taus. Since an 
event typically has 8 -fold solutions, we are free to keep all of the real solutions for an event and drop the ones that 
are complex. Note that we have assumed the input masses uncertainties for f and to be around 5% in event 
reconstruction. We find the result that with the inclusion of detector resolution effects that around ^2/3 events have 
only complex solutions, and we are left with ^ 30 stau pair events with one or more real solutions for the study of 
the new particle’s spin. 

Next we study the discriminating power of As after the detector resolution effects are included to establish whether 
it is still possible to distinguish the spin of the new particle spin using such a limited number of events. Simulating 
with a very large number of events including detector resolution effects we can deter mine that Aferm ion = —0.048 
and ^scalar = 0.160. The corresponding statistical uncertainty can be estimated as ^(1 — ^scalar)/-^- With 30 ^ 
1.2 X (1 — Agj.,j[g^j.)/(xlscalar “ ^fermion)^, we Can conclude that the ILC with integrated luminosity of 1000 fb ^ will 
only be able to distinguish the particle spin by more than Ict for our benchmark process . 


^ More than Scr can be achieved at the final state of ILC with integrated luminosity of 8 ab' -1 m 
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3. Tau polarization 

Since we are studying the polarization through r —>■ ttv channel and its decay branching ratio is ^ 11 %, we will 
have in total 


Np ~ lOfb X lOOOfb”^ X 0.72 X (11% x 2) - 1100 (56) 

tau decays for our impact parameter analysis. One of the advantages of the impact parameter analysis is that it does 
not suffer from any failed reconstruction problems. Furthermore, the ILC can have a quite good resolution for the 
impact parameter measurement. 

We hrst study the detector effects by using a very large number of signal events, since the statistical problem can 
be studied separately. The corresponding Aqip for different polarization of tau after considering the detector effects 
are given in Tab. |TTj 


Pr 

1.0 0.9 

0.2 0.3 0.1 

-1.0 -0.9 

Agip 

0.698 0.685 

0.605 0.615 0.591 

0.478 0.466 


TABLE II: Aqip for different degrees of tau polarization after detector resolution. 


We can conclude from the table that due to the excellent impact parameter resolution of the ILC, Agip works 
essentially as well as before even with detector resolution effects included. A 0.1 variation of Pr will typically lead 
to 0.01 changes in Aqip, by estimating its statistical uncertainty as y^(l — Aq^p)/N we can conclude that 1100 tau 
decays will be able to provide the measurement of tau polarization within ^ 0.25 precision. 

Note that we have assumed full purity of the tt channel for studying tau polarization. Experimentally the tt and 
p channel can be distinguished using the energy sum of charged and neutral particles. A dedicated study [S5] using 
neural networks has shown that the purity of the tt mode can reach 96%. As for the dcip distribution of the p channel, 
we find it is almost the same for all different tau polarizations with Aqip ~ 0.61, i.e., it is the Aqip corresponding to 
Pt ^ 0 in TT the channel. So, the small residual admixture of the p channel with the tt channel will act to reduce the 
estimation of the tau polarization by only a few percent. 


V. CONCLUSION AND DISCUSSION 

In this work we have exploited the relatively long lifetime of the tau lepton to demonstrate a new way to extract 
new physics parameters. 

For a tau undergoing a 3-prong decay, its direction can first be determined from the location of its displaced vertex 
and subsequently its energy can be determined up to a two-fold ambiguity. We saw that the distribution of the 
false tau energy solution of this channel is insensitive to the tau polarization. By extracting the end point energies 
(Amax/min) from the reconstructed tau energy distribution, we showed that the mass reconstruction precision of new 
particles (r and x in our case) can reach a precision of ~ 10 GeV with 3000 three-prong tau decays. 

If we require both tau leptons in the final state to go through 3-prong decay such that we know both tau directions, 
then the whole system can be reconstructed. Even though there is an eight-fold ambiguity, we find the distribution 
of those false solutions have a somewhat flatter shape than the true solution distribution, i.e., more concentrated in 
the central region for the scalar and more concentrated in the forward/backward region for the fermion. By studying 
the statistical uncertainty and using the distribution of all of the eight-folds solutions, we find that only a relatively 
small number of stau pair events (^ 120) are required to establish a 3tT differentiation between scalar and fermion 
final states. 

We have also proposed a new method to measure the tau polarization in the t ^ ttv channel, i.e., by using the 
impact parameter distribution of the charged pion in the final states. This method has the advantages of being 
more easily accessed experimentally and of being measurable with high precision at the ILC. We also find the impact 
parameter distribution is sensitive to the tau polarization while being very insensitive to the tau energy. This was 
seen to be particularly true in the parameter region of interest. We observed that the impact parameter distributions 
for different tau polarizations appear to intersect at approximately dcip 41/rm for our signal processes. 

With the assumption of a relatively pure right-hand polarized beam, we find that the backgrounds and NLO 
correction effects are typically more than one order of magnitude smaller than our leading order process and so these 
effects have been are neglected in this first analysis for simplicity. Assuming a signal production cross section of 10 
fb and taking into account realistic detector resolution effects we have found that with integrated luminosity of 1000 
fb“^ the mass reconstruction precision can reach ^ 20 GeV for our benchmark point {rrif = 300 GeV, = 50 
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GeV). In attempting spin reconstruction we encounter 3 quadratic equations. Taking into account the effects of finite 
detector resolution we find complex solutions in some cases and these must of course be rejected. The rejection rate 
for false solutions was found to be much higher than that for true solutions and this leads to a welcome increase in 
the ratio of true to false solutions. This then leads to an increased difference between the superposed distributions 
for all real solutions for both scalars and fermions. After taking into account all detector effects we conclude that 
200 reconstructed stau pair events would be enough to resolve the spin of a new particle to 3tT C.L. Since the ILC 
would provide a good resolution for the impact parameter, the discriminating power of doip after taking into account 
detecter resolution effects is essentially unchanged. The ILC with integrated luminosity of 1000 fb“^ is expected to 
be able to resolve the tau polarization to within a precision of ~ 0.25. 

A comparison with some previous studies |35[ 1391140j shows that this earlier work reported better resolution of 
the stau mass and tau polarization than the results presented above. However, this is primarily due to the limited 
statistics that result from our choice of our relatively heavy benchmark point, i.e., ruf = 300 GeV. For example, in 
Ref. [33] there were ~ 9000 stau pair events before any selections were used for f 2 and, in addition, it was assumed 
there that the neutralino mass m^o was already known. With these assumptions, the end point Flmax can be fitted 
with an uncertainty of ~ 2 GeV, which leads to an uncertainty in rrif^ of ~ 10 GeV. Moreover, as pointed out in 
this reference, the reconstructed rrif^ is very sensitive to the presumed value of e.g., an error of 80 MeV on m^o 
translates into an additional error of 1.4 GeV on rrif^- Our quoted lower resolution represents a conservative estimate 
without assuming a known value of m^o and also in part results from a smaller statistical sample in the 3-prong tau 
decay channel of our benchmark point (2600 three-prong tau decays before selection). As for the tau polarization 
measurement using the tt spectrum quoted in Ref. [33], they used 7.92 x 10^ stau pair events before selection for rifi 
with a corresponding uncertainty on the tau polarization of AP,- ~ 0.1, which reduce to 0.06 after considering the 
correlation between fitted normalization and the polarization. Assuming the same number of signal events in our 
study the uncertainty in Aqip is i5Agip ~ 0.0076. As can be seen from Table [III the corresponding change in the tau 
polarization is AP-r < 0.1. In summary, we see that a similar precision in the tau polarization measurement can be 
achieved with the same number of assumed events. 

In summary, this work has demonstrated an approach to searches for new physics particles that exploits the relatively 
long lifetime of the tau and the resulting displaced secondary vertex of the tau decay. We have shown that these 
techniques allow a determination of the mass and spin of the new physics particle. They also provide information 
about the couplings of the new physics particle that can be inferred from measurements of the tau polarization. 
This work represents a valuable complementary approach to the determination of these quantities and the precision 
obtained is comparable with other approaches. 
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